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Q\ , The eigenstates and the scattering transmission for two interacting elec- 
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trons are found exactly for / quantum dots, including the hybridization with 
the states on the leads. The results imply limitations on the validity of 
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the Coulomb blockade picture. The ground states for 1 = 1, 2 on a one 

O 



dimensional chain (modeling single and double quantum dots) exhibit quan- 
tum derealization and magnetic transitions. The effective transmission T 
of two interacting electrons through one impurity (/ = 1) is enhanced by 
a renormalization of the repulsive interaction, when one of the electrons is 
captured in a strongly localized state. 

I. INTRODUCTION 

There has been much recent interest in the effects of interactions on the localization of 
electrons in disordered systems [|TJ. Following many recent papers, we address this issue 
for the simplest case of two electrons. Recent numerical work 0-U found that interactions 
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help to delocalize the electrons, yielding a metal-insulator transition. Instead of discussing 
a fully random system, we consider a dilute system with only I (out of N) impurities. We 
then discuss the 'static' eigenstates and scattering situations. In Sec. 2 we show that the 
exact eigenenergies of the two electrons, in a tight binding model with on-site interactions, 
are found from the eigenvalues of a (small) I x I matrix, which involves the solutions of the 
le Hamiltonian. This yields a new calculational scheme, which should be very effective for 
numerical studies. 

Using this scheme, we obtain consequences for quantum dots f|, and our discussion 
emphasizes the importance of the hybridization of the wave functions with the 'band' states 
of the 'leads'. Quantum dots have potential applications as artificial atomic or molecular 
devices, and are very instrumental in studying strong correlations and their implications. 
Much of the discussion in the literature used the Coulomb blockade approach ||, which 
asigns a constant repulsive energy for each pair of electrons on the dot. In contrast, we 
find that the interacting eigenenergies are bounded between consecutive non-interacting 
two-electron (2e) energies, so that the energy cost due to the interactions is usually much 
smaller than the average on-site repulsion (U). 

In Sec. 3 we apply our scheme to specific models of quantum dots. A single quantum 
dot coupled to electrodes is modeled by one impurity (/ = 1) on a one dimensional (ID), 
iV-site chain 0. This represents a special case of the Anderson model ||, with a specific 
momentum-dependent hybridization between the impurity and the conductance band. We 
find that the behavior of the two electrons on the "dot" has a rich phase diagram, as 
function of the "dot" site energy eo and the hybridization. These parameters can be tuned 
experimentally, by varying the voltage on a gate coupled capacitively to the dot and the 
barriers between the dot and the leads ||. We also find a derealization transition at large U, 
for sufficiently negative eo, giving an analytic demonstration of the metal-insulator transition 
mentioned above. More complex examples include a cluster of / impurities coupled to ID 
leads, representing one large dot (JTTJ, two separate impurities, representing double quantum 
dots (or artificial "molecules") Ol, etc. 



Section 4 generalized the above discussion to scattering situations. For a simple model, 
with I = 1, we find a large enhancement of the transmission by the interactions, reflecting 
again derealization by interaction. We find peaks in the transmission (or in the conduc- 
tance) as function of the gate voltage on the dot. These peaks have nothing to do with those 
discussed in the Coulomb blockade literature. 

II. SCHEME FOR / IMPURITIES 

We start with the le tight binding Hamiltonian 

i 

7~Lo= (t nm \n)(m\ + h.c.) + ^2ei\i)(i\, (1) 

(n,m) i=l 

where |i) is a (spin-independent) state localized on site i, the first sum runs over all the 
site pairs in the system (including the impurities) and the second sum runs over the / 
impurities. We first find the le eigenstates \a) = J2 n (pa(n)\n) and eigenenergies e a of TCo. 
For the interactions, we consider only on-site repulsive interactions, TCi n t — ^f=i^i^i,T^»,l- 
This interaction has no effect on the triplet states, which vanish when both electrons are 
on the same impurity. Therefore, we discuss only singlet eigenstates, ^f(i,j) = with 
electrons at sites % and j. 

Rewriting the eigenvalue equation (7i — E)^f(i, j) = in the form \I/ = (E — Ho)' 1 ^^^ , 
with 7i = 7i + 7^int, we can now easily obtain the equations Aj = Hj =1 SijAj, where 
A; = ^yUi^(i } i) } and Sij/^JUiUj = G%{ii,jj) = (ii\(E — 1 1 ) is the non-interacting 2e 
Green's function in which the two electrons are on the same site, cf. ||. The new eigenvalues 
are thus found by requiring that the I x I determinant D(E) = 1 1 Sy — Sij\ \ vanishes. Finding 
the / eigenvalues Si(E) of the matrix S, the new (^-dependent) exact eigenenergies {E} 
solve the equations Si(E) = 1. The Aj's are then given by the eigenvectors related to Si(E), 
and \l/(m, n) = HiG^mn^ijy/TTiAi. 

Using the non-interacting singlet basis, ^Q b (i,j) = [0 a (i)0 b (j) + <fi a (j)(f) b (i)]/2 ( - 1+Sab ^ 2 , 
with eigenenergies e a b = e a + e?, and e a < e&, one has 
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G%{ijM) =Y,Ui)U3W a {k)m)/{E - e ab ). (2) 

ab 

Since each Sij has poles at every e ab , each Si will generally also have such poles. Ignor- 
ing special symmetric cases, when the residues of some of these poles vanish, one has 
oc 1/(E — e ab ) when E « e ab . In this approximation, D{E) « (— 1) 7 (1 — TrS) ~ 
(_!)/[! _ (U) ab /(E - e ab )}, with (*7) a6 = E^I^MI 2 - As £ crosses e ab , (-1) I+1 D(E) 
jumps from — oo to oo. As E increases between two consecutive non-interacting eigenvalues, 
(— 1) I+1 D(E) varies smoothly from oo to — oo. Thus, D(E) = must have at least one (and 
up to /) solution(s) between every pair of such consecutive energies, and the new energies 
maintain the sequence of the non-interacting ones. For example, the cost of adding two 
electrons into the ground state, A gg = E gg — e gg , is thus bounded by the le distance to the 
first excited level, e u — e g , which may be much smaller than the Coulomb blockade value 
(U)\ It is interesting to note that small (and not evenly spaced) values of A were observed 
in a series of 2D quantum dots [12|], apparently contradicting the simple Coulomb blockade 
picture. It is tempting to relate these observations to our result. 



III. ONE AND TWO IMPURITIES: GROUND STATE AND PHASE 

TRANSITIONS 

The simplest example concerns one impurity (or "dot") on a ID wire, closed with periodic 
boundary conditions. We place the impurity at site % = 0, with energy eo and with real matrix 
elements t 01 = t at_i = t . The other nearest neighbor matrix elements are set at t n — t, 
for n = 1,2, ...,N — 2. The le spectrum contains N/2 odd eigenstates, <fik{n) oc sin(fcn), 
which are not affected by the impurity, and N/2 even ones, whose energies are shifted. These 
even solutions exhibit a rich phase diagram in terms of e (in units of t) and 7 = (to/t) 2 , 
as shown in Fig. [1] |TJ|]: In region A all the states are delocalized, with almost unshifted 
energies, in the 'conduction band'. In region C (or D+F) there exists one bound state above 
(or below) the conduction band. Finally, both bound states exist in region B. An important 
detail concerns the normalization: except for the trivial case when 7 = 1 and eo = 0, the 
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amplitude of the 'band' state |0fc(O) | 2 is found to be proportional to sin 2 k, vanishing at the 
band edges. This will turn out to be crucial below. 

For 1 = 1, adding the interaction U\ = U, D(E) = Sn(E) — 1. Indeed, Sn(E) jumps 
from — oo to oo as E crosses each non-interacting energy e a b, and we find one new eigenvalue 
between every pair of non-interacting energies, as described above. In regions B, D and F, 
the non-interacting ground state has two bound electrons, with energy e gg = —2\e g \, and 
the new ground state energy is always bound between this energy and that of one bound 
and one 'free' electron, e gu = —\e g \ — 2t. For finite N, Sn would diverge towards — oo as 
e gu is approached, implying a persistent doubly bound solution below the band for all U. 
However, when iV — > oo we replace the sums over the band states by integrals. Since 10^(0) | 2 
vanishes at e gu , Sn(t g u) = Us c remains finite. In region D one has s c < 0, so that D(E) = 
has a discrete doubly bound state for all U. However, in region F one has s c > 0, implying 
a disappearance of this bound state for U > U c = 1 / s c ; the ground state energy is now at 
the bottom of the band, implying an 'ionization' of one electron and an 'insulator to metal' 
transition from D to F. This transition occurs at a lower value of U for larger negative 
values of eo, when the localization length of the bound electrons is smaller, and they feel the 
interaction more strongly. 

It should be noted that although the bound state is a singlet, with total spin zero, the 
new ground state in region F has one bound electron and one "free" electron. Such a state 
does not feel the e-e repulsion, and is thus practically degenerate with the slightly lower 
triplet state (for large N, the difference is of order 1/N). Unlike the "insulator" singlet (or 
"antiferromagnetic" ) ground state, which has no net magnetic moment, this "metallic" state 
in region F is paramagnetic. This difference should be measurable in an external magnetic 
field. Our model also exhibits interesting finite size effects for N < oo, and interesting 
antibound 2e excitations above the band [fn}| . 

The analysis is somewhat more complicated for 1 = 2, where the 2x2 matrix S has 
symmetric and antisymmetric eigenvalues, corresponding to products of even-even (and 
odd-odd) and of even-odd le states, respectively. The resulting phase diagram now also 
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depends on the inter-impurity distance R, and the region equivalent to F in Fig. [I] shrinks 
to a small 'bubble' 113 . 



IV. ONE IMPURITY: TRANSMISSION 

So far we emphasized only the ground state, for a problem with periodic boundary 
conditions. In real 'quantum dot' experiments one is more interested in scattering situations, 
when electrons are sent from one side and one calculates the current or the transmission 
through the dot. For the simplest exact calculation of the transmission for two interacting 
electrons, we replace the above tight binding Hamiltonian by its continuum limit, 

H(xt, ai, x 2 , o- 2 ) = Ho(x 1 ,x 2 ) + U8(x 1 )5(x 2 )S ai - a2 , 

Ho(xi, x 2 ) = Hspixt) + H sp (x 2 ), (3) 

where X{ and crj are the coordinate and the spin component of the ith electron, and 7i sp is 
the single-particle Hamiltonian, independent of the spin components. 

Again, we consider only the singlet spatially symmetric wave functions, $f(xi,x 2 ) = 
^(x 2 , xi). At total energy E, we split \I> into ^(x%, , x 2 ) = ^0(^1? x 2 )+^s{xi, x 2 ), where ^0 is 
the solution of Ho(xi, x 2 ), with the same energy E + = (E+iw) (with r\ — > 0), {Ho — E + )^/ = 
0. For the on-site interaction (^|) it then follows that ^/s(xi,x 2 ) = UGe(xi, x 2 ; 0, 0)\l/o(0, 0), 
where Ge is the two-particle Green's function of the interacting Hamiltonian. 

For the model Hamiltonian given by (||) one can express the Green's function of the 
interacting system, Ge, in terms of G° E , discussed in the previous section. The definitions 
of these two Green's functions yield the Bethe-Salpeter equation 

G E (x 1 ,x 2 ;x' 1 ,x' 2 ) = G%{xi 1 x 2 ;x' l) x' 2 ) + UG° E (0,0;x' 1 ,x' 2 )G E (x 1 ,x 2 ; 0,0), (4) 

and hence G E {x 1: x 2 ; 0, 0) = G E {x u x 2 , 0, 0)/[l - UG° E (0, 0, 0, 0)]. Thus, ^{x u x 2 ) = 
Mxi,x 2 ) + F B Gk(xi,x 2 ,0,0)tf (0,0), where F E = U/[l - £76^(0,0,0,0)]. Thus, V s is 
determined solely by the non-interacting Hamiltonian 7io- 
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We next calculate the quantum average of the current density operator at x = xo, in the 
exact singlet state 

2eh r d 
j(x ) = J dx 1 dx 2 5(x 1 - x )V*(xi,X2)—i$?(xi,X2). (5) 

The explicit calculation of j now requires only integrals involving the non-interacting func- 
tions ^0(^1^2) and non-interacting states <p p (x). In what follows, we shall assume that 
^0 is given by the singlet combination \I>o 9 , as defined in the previous section, and that the 
total energy is given by E = e p + e q . 

To proceed, we choose a simple 5-function attractive potential, 

which has one bound state, <pb = y/ne~ K ^ with the inverse localization length k = mV/h 2 , 
with eigenenergy — e& = — h 2 K 2 /(2m), and "band" scattering wavefunctions <p p = (e ipx + 
r p e ip ^)/\^L, with the reflection and transmission amplitudes r p = in/ (p — in), t p = p/(p — 
in), and with eigenenergy e p = h 2 p 2 / (2m). 

There are two physical situations which are of interest for the non-interacting wave 
function ^/ p ) q . The first corresponds to two propagating electrons, impinging from the left, 
when both p = p 1 and q = p 2 represent "band" states. In this case, the current density of a 
'macroscopic' system (that is, for L — > 00) is unaffected by the interaction U, and remains 
as in the non-interacting system, j = eh(pi\t pi \ 2 + p 2 \t P2 \ 2 ) / (mL) + 0(U/L 2 ). 

The second, more interesting, choice for ^q 9 arises when one electron is propagating, with 
p representing its wave vector, and the other is captured in the bound state, i. e. q — b. Now 
the total energy is E(p, b) = e p —e b > —e b . In this case, the terms coming from the interaction 
are of the same order as the non-interacting ones. A long but straightforward calculation now 
yields j = ehpT (p) / (mL) , where the effective transmission is T(p) = \t p \ 2 — 2m^R.(F E r p t p ) /h 2 . 
Remembering that F E = U/[1- UG° E (0, 0; 0, 0)], the new term in T will have a 'resonance' 
when KF^ 1 = U^ 1 — ?R,G E (0, 0; 0, 0) = 0. This equation is similar to the equation Sn = 1 
encountered in the previous section. 



For an explicit calculation of G° E (0, 0; 0, 0), we introduce an upper cutoff W = h 2 u! 2 /2m 
on the "band" states. We then found that dtG E (0, 0; 0,0) < for small k, and 
$IG E (0, 0; 0, 0) > for large k. In the former region there is no 'resonance'. In the lat- 
ter region there can be a resonance for sufficiently large U. This yields a very interesting 
dependence of the effective transmission T on the single-electron parameters and on U. 
Specifically, Fig. |2] shows the dependence of T(p) on k, for a finite value of U and for two 
values of p. Clearly, T(p) increases significantly as the result of the interaction, reflecting 
derealization due to the interaction. At k = one has T — 1, independent of U. At small 
K, T always decreases, due to scattering by the single electron potential, following the non- 
interacting case. However, U effectively screens this attractive potential, and at small p this 
results in a double peak structure of T, as shown in the figure. Qualitatively, these peaks 
arise when the screening cancels the attractive potential. At larger p the minimum between 
the peaks grows, the two peaks join into a relatively broad plateau which decays at large p. 

This project is supported by the Israel Science Foundation and by a joint grant from the 
Israeli Ministry of Science and the French Ministry of Research and Technology. O. E-W 
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FIGURES 




FIG. 1. eo — 7 phase diagram for the single impurity case. 
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FIG. 2. The transmission T as function of k (in units of the cutoff u), for U = 5h 2 /(2m). The 
solid line is T(p), the dashed line shows \t p \ 2 . (a) p = 0.04a;, (b) p = O.loo. 
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